Phase separation in the large-spin t- J model 

Antimo Angelucci and Sandro Sorella 
International School for Advanced Studies, via Beirut 4, 34014 Trieste, Italy 



We investigate the phase diagram of the two dimensional t-J model using a re- 
cently developed technique that allows to solve the mean-field model hamiltonian 
with a variational calculation. The accuracy of our estimate is controlled by means 
of a small parameter 1/q, analogous to the inverse spin magnitude 1/s employed in 
studying quantum spin systems. The mathematical aspects of the method and its 
connection with other large-spin approaches are discussed in details. In the large-g 
limit the problem of strongly correlated electron systems turns in the minimization of 
a total energy functional. We have performed numerically this optimization problem 
on a finite but large L x L lattice. For a single hole the static small-polaron solution 
is stable unless for small values of J, where polarons of increasing sizes have lower 
energy. At finite doping we recover phase separation above a critical J and for any 
electron density, showing that the Emery et al. picture represents the semiclassical 
behaviour of the t- J model. Quantum fluctuations are expected to be very important 
especially in the small J - small doping region, where phase separation may also be 
suppressed. 
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I. INTRODUCTION 



The discovery of high-temperature superconductors renewed considerable interest and 
attention to the study of strongly interacting electron systems in one and two dimensions, 
because it is widely believed that the anomalous properties of such materials may be related 
to the strong Coulomb repulsion and to the low effective dimensionality. 

One of the most interesting models which exactly incorporates the constraint of strong 
Coulomb repulsion is the well known t- J model, which will be the main topic of the present 
paper. Several approximate techniques have been proposed so far to deal with this model 
hamiltonian. We mention, for example, the self consistent Born approximation [l], the semi- 
classical approach 0, and recently the limit of infinite dimensionality ||, which have been 
developed for the one-hole case. In the general case of arbitrary density, the large- N expan- 
sion with slave-boson and the slave-fermion techniques j|] should also be mentioned, as well 
as the vast amount of numerical and variational works ||. The results of all these different 
methods are quite controversial, especially concerning the question about the presence in 
the model of superconductivity and/or of marginal Fermi liquid behaviour. 

Quite recently Emery et al. tried to explain the full phase diagram of the model, 
basing their analysis both on a variational argument in the small J region and on some 
numerical evidences emerging from a 4x4 lattice exact diagonalization. They speculated that 
phase separation should occur throughout the full phase diagram: the electron-system phase 
separates into a phase where the holes move in a fully polarized state and into an electron-rich 
phase characterized by a Neel antiferromagnetic spin-order. Accordingly, phase separation 
occurs below some critical doping <5 C , where 5 C — > for J — ► 0. Castellani et al. then argued 
that superconductivity may occur close to the phase separation boundary. A numerical 
work by Ogata et al. || provided the complete determination of the one dimensional phase 
diagram and evidenced this property in one dimension. 

After the Emery's work, Putikka et al. performed a systematic high-temperature 
expansion on the t- J model and found evidence that the separated phase should appear only 
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for large values of J. Although this work is surely non-conclusive, they pointed out that 
there is not convincing evidence that the separated phase should be continuously connected 
through the whole phase diagram, because the small J and the large J region should behave 
in a very different way. 

Due to the present controversy, we derive here a consistent mean-field calculation on the 
t- J model using a recently developed technique [[H| allowing to control the accuracy of the 
mean-field estimate by means of a small parameter 1/q, which is similar to the parameter 1/s 
- the inverse spin magnitude - used in the spin-wave theory of quantum spin systems. This 
technique has the following advantages: i) does conserve the symmetries of the t-J model 
hamiltonian, ii) takes exactly into account the constraint of no-double occupancy, and in) 
ensures that the classical estimate of the energy is variational for any value of q. In this 
way we generalize a very useful property known to hold for spin systems, i.e., the classical 
solution is independent of the spin magnitude and the corresponding energy is therefore 
variational. 

With our mean-field approach we are able to reproduce the phase separation over the 
small and the large J region, supporting at the semiclassical level the Emery et al. picture. 
We have made an intensive and systematic numerical work in order to calculate the true 
phase diagram at the mean-field level without imposing any a priori order parameter. At 
the end the phase diagram looks extremely simple, but unfortunately poor. However our 
calculation represents only the first step towards a complete determination of the phase 
diagram. As standard in any semiclassical approach, the second step would be to include 
quantum fluctuations about the mean-field. With the present technique corrections to the 
classical energy can be introduced systematically. 

The paper is organized as follows: in Sec. |I| we review the graded Holstein-Primakoff 



map [|1(J for the graded algebra spl(2,l) and present the associated coherent states. These 
mathematical tools are developed starting from the observation that spl(2,l) is the algebra 
of the operators entering the strongly correlated electron systems and that it has a series of 
representation characterized by a parameter q analogous to the spin magnitude. The small 



parameter 1/q allows to systematically define distinct large-spin limits of the t- J model. 
In Sec. [TTT] we investigate the physical implications underlying the different generalizations 
of the model and show that previously proposed effective hamiltonians, notably the Kane- 



Lee- Read and the spinless fermion hopping hamiltonian |TI| , can be derived as particular 
cases of our approach. In Sec. [TV] we then derive the variational total energy functional which 
solves our mean-field theory and present the results concerning the numerical investigation. 
Sec. [V] contains the conclusions. 

II. SPL(2,1) COHERENT STATES 



The t- J model hamiltonian [12| is defined by 

H tJ = tJ2 ^iAc aj + J -x E ^' (SiSj - ijvA) , (2.1) 

ij,a ij 

where a = 1, 2 are spin up and spin down indices, respectively; fiy is the adherence matrix 
with periodic boundary conditions connecting nearest neighbors sites, i, j denote the lattice 
coordinates, and 

Cli = (1 — C*2i c 2i)Cli , c\i = cjj(l — C2jC 2 j) , 

C 2i = (1 - c\ i C X i)c 2 i , 4 = - c li c li) , ( 2 - 2a ) 

are hole creation/annihilation operators and 

= (cJjCii + C^i c 2i) , S 3i = -(cJjCij — cliC2i) , 



Sli — ~{Ai c 2i + cliCu) , S 2 i — — {c\iC2i — cliCu) , 



(2.2b) 



are charge and spin operators, where c| i; c' 2i are the electron operators, and doubly occupied 
sites are excluded. The exchange constant is always positive: J > 0. 

The 3x3 Hubbard matrices are defined by taking the expectation values of the operators 
(2.2) in the restricted single-site Hilbert space |0j), = c|j|0j), |2j) = c^jOj) 



Xai = (Oi| C ai , d = (Oi| iVj | A) , 

X?=<Oi|4lA) , I A) , a,/3 = 0,l,2 



(2.3) 



The matrices an d Xai are conjugate - i.e., Xi — Xai- However for later convenience we 
prefer to use the notation with the raised index. The matrices ( p. 3D are useful because the 
constraint of no doubly occupancy can be automatically taken into account by replacing in 
Eq. ( |2.1|) the operators (2.2) with the corresponding matrices Q2.3Q . Defining Q 0i = (J— |Cj), 
where I is the identity matrix, the Hubbard matrices (|2.3j) can be divided into odd (i.e., 
"fermionic") generators, Xai, Xii an d even (i.e., "bosonic") generators Q^i = (Qoi,Qi) - 
which we collect in four-vectors - satisfying the commutation/anticommutation relations of 
the spl(2,l) graded algebra flT0| , |T3"l 



[QiihQvj] — dijitOfiv Q\j , (2.4a) 
\xS,Qh] = 8vI(*M, (2.4b) 

X"} = 8iA°*)lQri > iXai, Xbj} = , (2.4c) 



where the completely antisymmetric tensor is normalized by eoi2 3 = 1; whereas the spin 
four- vector is built from the standard Pauli matrices as (cr^)* = (5 a b,<?ab)- The greek four- 
vector indices are raised and lowered using the metric tensor g^ v = diag(l, —1, —1, —1), and 
a summation over repeated "Lorentz" and latin ( "spinorial" ) indices is understood. The 
generator Q 0i has to be introduced in order to close the spl(2,l) algebraic rules ( |2.4b|J2.4cD . 
For the time being our considerations will refer mostly to single-site quantities and to simplify 
the notation the site index will be dropped - if not otherwise needed. 



As explained in details in Ref. ||I4|| , the several classes of irreducible representations of 
the spl(2,l) graded algebra are labelled by the eigenvalues of the operators Q , Q 2 , and Q 3 , 
respectively denoted qo, q(q + 1) (q is called isospin, and is an integer or half-integer number), 
and q%. The basis vectors are thus denoted with \qo,q,q3). Among them, those relevant for 
our investigation are the so called atypical representations. They are characterized by a linear 



relation between q and q and are thus identified only by the value of the isospin q - this 
being the parameter which we shall use analogously to the spin magnitude s in the standard 
spin-wave theory of the Heisenberg hamiltonian. In each such atypical representation, of 
dimensionality (4g + 1), the basis vectors can be grouped into two multiplets of (2q + 1) and 
2q states with quantum numbers 

\Q,Q,Q3), ?3 = ~q, -q+ l,---,q, (2.5a) 
11. 1 1 . , . 

\q+ 2' 9 - 2' g3 '' Qs = ~ q+ 2''"' g ~ 2' ^ ' 

which we assume as even and odd, respectively (the grading of the states is a pure matter 
of convention). The relative normalization of the two multiplets is not a priori fixed. For 
the isospin value q—\ the states (2.5) are the real spin-^ and hole (spin-0) states, respec- 
tively, and the matrices defining the representation turn in the 3x3 Hubbard matrices. 
Thus the operators (2.2) entering the t- J model belong to the fundamental q = \ atypical 
representation of spl(2,l). According to Eqs. (2.5), models of strongly correlated electron 
systems can be generalized in a natural way by enlarging the dimensionality of the spl(2,l) 
representation. In particular we are naturally lead to identify the multiplet ( |2.5a|) as the 
generalized spin state, and the multiplet ( [2.5b]) as the generalized hole state. The different 
value of go between the two multiplets is related to the grading. Defining the "fermion 
number operator" F = 2Qq — 2ql - where / is the unit matrix of the g-representation - 
then the spin multiplet ( |2.5aj ) has F = and the hole multiplet ( |2.5b| ) has F — 1. Thus for 
a many-particle system the conservation of the number of particles can be rewritten as the 
conservation of the fermion number 

N N 

Y.^ = £(% - 2gJ) = N h = N- N el ; (2.6) 
i=i i=i 

where N e \ is the number of electrons, iV is the number of sites, and thus is the number 
of holes. 

Assuming the maximal isospin state as the reference vacuum, |0) = \q,q,q) (not to 
be confused with the electronic empty state), and introducing a canonical boson a and a 
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canonical spinless fermion c, then the (4g + l)-dimensional Hilbert space spanned by the 
states (2.5) can be put into correspondence with the Fock space generated by the states 
\B, n) = (a^) ra |0) [corresponding to Eq. ( |2.5a| )], and \F, n) = c^(a^) n \0) [corresponding to Eq. 
( [2 -5 b| ), with \F,0) = \q + ~, q — ~, q — ^)], satisfying the operatorial relation a)a + c^c < 2q. 



In this basis, the generators (2.4) can be represented as follows pUfl 

Q = q + ^c, xi = c^2q - a)a - dc , 

Q 3 = q — a^a — -c'c , x 1 = \j2q — a^a — etc c , 



Q + = J2q — a^a — tfc a , \2 = c T a 



Q_ = <Jyj2q - ata - etc , % 2 = a f c , (2.7) 
where Q± = Qi ± «Q2- Note that the representation of Xi (and x 1 ) can be simplified - i.e., 



we have xi — c'y 2<? — ata. The symmetric notation ( |2.7[ ) is useful to check that Eq. ( |2.7| ) 
satisfy to Eqs. (2.4). These equations represent the generalization to spl(2,l) of the usual 



Holstein-Primakoff (HP) transformation |L5] for su(2) and in this respect are similar to those 



of the Swinger bosons/slave-fermion representation [ IS] . Using the realization ( |2.7| ) of the 
generator Qo, the conservation law Q2.6| ) becomes the conservation of the spinless fermion 
number operator: Nh = J2i c l c i = Nh, so that at half-filling (i.e., N h = 0) Eq. (|2.7j ) reduce 
to the standard su(2) HP transformation 

Qs = s — a) a , 

Q+ = ^2s - ata a (2.8) 



Q_ = at y 2s — ata , 

where ata < 2s (using the standard symbol s instead of q), in agreement with the result 
that at half-filling the t- J model becomes the antiferromagnetic (AF) Heisenberg model. 

As for spin systems, we wish to describe spin configurations classically. To this aim we 
need the transformations (Q^,Xa) (P^, A a ) of the generators which preserve the spl(2,l) 
algebra. This can be conveniently done by choosing the original fermionic representation 
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(2.2) of the q = \ generators (2.3), and then considering the most general transformation 
of the electron operators which preserves the canonical commutation relations (2.4) as well 
as the constraint of no double occupancy. Simple algebra shows that the odd generators 
transform like su(2) operators in the fundamental representation: 

v i~ -i± ( it .9 _ii \ 
X x = e ll e l i I cos -e J 2 xi - sm -e 2 X2 I , 

■± ( 9 -± 9 _-± \ 

X 2 = e ll e l * I sin -e l 2 xi + cos -e '^a I , 

X l = X{, X 2 = X$, (2.9) 
whereas the even generators transform as 

Po = Qo , Pk = RkmQm, (2-10) 

where Rkm is the SO (3) rotation matrix, whose explicit realization can be easily derived by 
means of the anticommutation rule ( ^.4c ). By retaining only the leading order terms of the 



HP transformation (|2.71 ) we have Q = Q 3 — > q, Q± — > 0, xi ~ ¥ \/2<Z c *> X2 —> 0, so that in 



the large-g limit the rotated even generators (|2.10|) become classical c-numbers 

P[i -> 9(1, n) =p^ , 
n = (sin 9 cos </>, sin 9 sin 0, cos 9), (2- 11) 

and the rotated odd generators (|2.9|) become proportional to the spinless fermion operator 



Xi -> fiqe^e-ii cos(^)c f , X 1 = X\, 

X 2 -> fiqe^e^ and) J , X 2 = X ] 2 , (2.12) 

where we have set u> = 7 + ip/2. 

In the previous equations 9 and are the physical angles which parametrize the classical 
spin configuration [see Eq. ( [2. lip ]. Instead, the third Euler angle ip and the global fermionic 



phase 7 are unessential. In fact, ip can always be eliminated by redefining Q± = e Tl ^Q±, 
X2 = e~ l ^/ 2 X2i Xi = e 2 Xi> an d 7 can be arbitrarily fixed by changing the relative normal- 
ization of the two multiplets (2.5). 
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Similarly to spin models, we can give to Eqs. ( |2.11| , pT2j ) a quantum mechanical meaning 



by introducing the even coherent state of the generalized spin multiplet ( |2.5a| ) and the 
odd coherent state \ flf) of the generalized hole multiplet ( |2.5b| ) 

\Q b ) = (cos^) 2 *exp(tg~e^Q_)|£,0) , (2.13a) 
lfy> = (oos^^eapCtg^g^lF.O) . (2.13b) 



The states (2.13) are orthonormal and by generalizing standard results [[IT]] on su(2) coherent 
states we find 

(n b \Q li \n b ) = q(i,ft)= Pli , (n b \Q li \n f ) = o, (2.i4a) 

(SlflQ^Slf) = (q+l } (q- l -)n\ = pi (2.14b) 



As for the su(2) coherent states, the classical c-numbers ( [2. lip are equal to the expectation 



value of the even generators between the even coherent states. The spin expectation value 
(Qf\Q\Qf) is in general finite and it is zero only for the value q — |, an obvious consequence 
of the fact that the generalized hole is actually a spin-(g — |) multiplet. Concerning the odd 
generators, it is easily checked that 

^-ft cos A |n/>, 



xiW = 




X^f) = 


^2qe 


X2\Mb) = 




x 2 \n f ) = 


\fzqe 




o, 



2 



v 2 /> - , 

^sinC^n/), 
-^+f) S in(~)|n 6 >, 

C a |n 6 ) = (2.15) 

where we have fixed the relative normalization of the states (2.5) by using the phase con- 
vention xi\B,0) = \/2qe % ^ u ~^\F,0). From Eq. ( |2.15| ) follows that the combination of odd 



generators 



X f = ( e if cos(-)xi + e _< » sin(-)x 2 J (2.16) 



creates the odd coherent state out of the even one, whereas the conjugate operator annihilates 
it 

X i\n b ) = \n f ), x\tt b } = (2.17a) 

We also have 

x\n f ) = \n b ), x t |fi / ) = o. (2.17b) 

The relations (|2.17a|) show that x an d X are the g-generalized hole creation and an- 
nihilation operators, respectively, and Eq. ( f2.17b| ) acquires the meaning of the parti- 



cle/hole transformation. Note that x > X are n °t spinless fermion operators, because for 
a generic g-representation their anticommutator (restoring the notation with the site index) 
is {xhxj} — $ij^(Qoj + Qjfij). However, this means that the identity 

(n|{ x I,xj}|n> = «5« (2.18) 

is satisfied, so that xti Xi ac t like spinless fermion operators when bracketed between any 
tensor product of even and/or odd coherent states |f2) = ® \£l a i) ( a = &>/)• 



III. LARGE-SPIN LIMIT 

Here and in the following we shall refer to the large-spin limit whenever the single-site 
Hilbert space of a strongly correlated electron system is enlarged by assuming a spin-g 
multiplet as the spin state and a spin-(g — |) multiplet as the hole state, and then one lets 
q — > oo. To develop the coherent states for spl(2,l) we have constantly used the operator 
Qoi, because this is the generator by which the spl(2,l) algebraic rules (2.4) can be closed and 
thus the spl(2,l) representation theory (2.5) applied. However in the t- J model hamiltonian 
there enters the charge operator iVj, and we have to define its large-g generalization. As a 
matter of fact, here we are faced with the more general problem of how to generalize the 
t-J model itself, because the large-g generalization of the charge and spin operators has a 
certain amount of arbitrariness. On one hand, to perform the large-spin limit of the t-J 
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model we may straightforwardly use the linear generalization of the operators (2.2). This 
amount to replace in Eq. ( |2~T| ) 

Cai^Xai, Si^Qi, Ni -> d = 2(2ql - Qoi), (3.1) 

where (xai, Q^%) have now arbitrary isospin magnitude q. The last modification generalizes 
the relation Cj = 2(1 — Q 0i ) of the original q = ~ theory. We have C^fi^) = 2g|fi bi ), and 
Cj|f2jj) = (2g — l)\Qfi), which means that the spin and hole states (2.5) have generalized 
spin - differing by half- unit - as well as generalized electric charge - differing by one unit. 

On the other hand, for the q = \ fundamental representation, and only in this case, the 
3x3 Hubbard matrices ( |2.3|) satisfy 



(3.2) 



Ci = (xlxu + x 2 iX2i), Qsi = -Axlxii - xlx2i), 



Qii = \ix\x2i + XiXu), Q2i = ^{x\x2i - xlxu), 



a consequence of the fact that whether the constraint of no doubly occupancy is enforced, 
the charge and spin operators ( |2.2b| ) are unchanged if we substitute the electron operators 



with the corresponding hole operators fl2.2a| ). According to the previous equation, we may 



think of the charge and spin operators (|2.2b ) as composite operators of the odd generators 



for any value of q, so that we may also define the large-spin limit of Eq. ( |2.1|) by means of 
the replacement 

Cai -> Xai, Si^ Qi = ^Xi^abXbi, Ni -> & = XiXai, (3.3) 

where Xai, Xt have now arbitrary isospin magnitude q. The present generalization amounts 
to consider Eq. ( |2.1| ) as a model of interacting holes and to preserve this dynamics for any 
value of the hole isospin. For any q ^ \ Eqs. ( |3.1| , |3.3| ) are not equivalent and in our analysis 
we shall investigate the physical implications underlying the two choices. 

The simplest generalization is obtained by choosing the linear realization ( |3.1| ). Using 
the arbitrarily rotated frame (P M , X a ) ( |2.1(j| , |2.9D , upon the substitution ( |3.1| ) the generalized 
t- J model hamiltonian becomes 
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ij,a ij 

= t JSJ QijXi X a j — — Qijg^ 1 Pfti.Pvj + 

ij,a ij,[iu 

MJN h , (3.4) 



where in the last identity we have applied the conservation law ( [2.6| ). Similarly to spin 
systems, using the linear generalization we preserve all the symmetries of the original model. 
Apart from a trivial constant, Eq. (|3.4j) is bilinear in the spl(2,l) generators, so that J = 2t 



remains the "supersymmetric point" |TJ| for any value of the isospin q. Replacing the HP 
map (|2.7| ) in Eq. (|3.4j ), in the q — > oo limit the potential energy becomes dominant, so that 
the large-spin limit of Eq. ( |3.4|) is given by the classical Heisenberg model 

H ( t f - = 3 -f E n« {uSi - 1) • (3.5) 

Z ij 

Thus for any electron density p = N e \/N the mean- field solution is always the Need spin 
configuration. We can easily understand this result noting that in the large-g limit the 
expectation values ( |2.14a|j2Tl4b| ) are equal 



lim = lim = fl,nj). (3.6) 

As a consequence in the large-g limit we still have two distinct states \ flu) and |fi/i) differing 
by the fermion number F^, however their spin and charge become equal; hence the mean-field 
theory turns in the AF spin hamiltonian, independently from the density. Thus in the large- 
spin limit we preserve all the symmetries but we miss the following important property of the 
exact model (|2.1| ): if the hole is localized on site io, the four bonds Jfij j do not contribute 
to the energy (we shall refer to this property as the "dynamics of the missing bonds"). 
This is a peculiar behaviour of the present large-spin limit. However all the developments 
are mathematically well definite and the expansion in fluctuations about the AF mean-field 
solution can be performed straightforwardly. We divide the lattice into two sublattices and 
in Eqs. ( [2.9| , |2.10| ) set 9 = = ip = 7 = on one sublattice, and 6 = ip = n , = 0, 7 = — 7r/2 
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on the other one. Then the replacement in Eq. ( |3.4|) of the HP realization (|2.7|) expanded 
in powers of 1/q gives the effective hamiltonian 
a T 

H KhR = £„ + y^ Qij (a\cLi + a]aj + a\a) + a^J 

ij 

-j2qtY,V lJ c\c j {a] + a i ) + ... (3.7) 
to describe fluctuations, where 

E = -ANJq 2 + 4Jq ^ cJq (3.8) 

i 

is the classical energy of the pure antiferromagnet plus a quantum contribution which shifts 
the chemical potential of the spinless fermions. In Eq. ( |3.7| ) we have displayed only the first 
few terms of the systematic expansion. They give exactly the hamiltonian proposed and 
applied by Kane, Lee, and Read (KLR) to describe the propagation of a single hole in a 
quantum antiferromagnet. In their approach the bosonic term and the three-body term - the 
latter responsible for the hole propagation through spin-wave emission and absorbtion - are 
treated on equal footing. This assumption corresponds to scaling t by a factor ^/2q and then 
define the large-spin limit keeping fixed the renormalized hopping parameter t = t/^/2q. 
Eq. ( p.7|) is well definite for any density and from Eq. Q3.8|) we see that in the large-spin 
limit ( |3.1| ) the dynamics of the missing bonds is present only when quantum fluctuations 
are taken into account. In particular setting q = |, we have E = —J{2N — AN h )/2, which 
is the energy of the classical antiferromagnet with delocalized holes. 

The dynamics of the missing bonds is the essential feature of the t- J model to understand 
phase separation in cuprate superconductors, and it is completely not accounted for by 
the large-spin limit ( |3.1|) . This property can be preserved if we resort to the nonlinear 
generalization (|3.3|) . All our considerations are simplified by introducing the projector onto 
the even sector 

^ = (2g+l)J-2Q oi , r^w) = |J2 W ), T i \n ft } = 0. (3.9) 

Employing the HP transformation ( j2.7[ ), it is easily shown that the even operators (Cj, Qi) 
of the set (|3.3| ) are related to the operators (C*, Qi) in Eq. ( [3.1|) via the projector Ti 
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d^rlQTi, Q i = TtQ i T i . (3.10) 

We collect the generalized charge and spin operators ( |3.3| ) in the four- vector = rfQ^Tj 
(note that we have the relation Cj = rfC^ = 2Q 0i ). Applying Eqs. ( |3~I||3.10| ) we then have 

(^u\Qfii\^u) =Pti,i, (^fi\Qni\^fi) = 0. (3.11) 

An independent check of Eq. (|3.11| ) can be obtained by using the definition ( |3.3|) and the 
rules Q2.15Q . Since the projector ( |3.9| ) commutes with the even generators Q^i it is easily 
proved that the su(2)xu(l) even subalgebra is preserved 

[Q/ui, Quj] = i5ije ^ x Q\j, (3.12a) 

so that the operators ( |3.10|) under the rotation (|2.9|) of the odd generators, Xai —> X ai , are 
(PftijXai), where 

Voi = Qoi = 2^°*' = X-iO k ah Xu = RhmQ.mii (3.12b) 



Rkm being the rotation matrix in Eq. fl2.10|) . Upon the substitution ( |3.3|) the generalized 
t- J model hamiltonian in the rotated frame fl3.12b|) eventually becomes 

ij,a 

^E^f^-7%V (3.13) 



2 2q ^ " J V J 4 

Here we have defined the model keeping fixed the rescaled exchange constant Jo = 2gJ, to 
weight on equal footing the kinetic and potential energy contributions (we explicitly display 
the dependence of H$-p] on the operators for later convenience). 

For arbitrary values of q we loose supersymmetry at the J = 2t point because the set 
of operators (Q^jXai) forms a representation of spl(2,l) only for q — |. However due to 
Eqs. (3.12) the model ( |3.13| ) remains rotationally invariant and, thanks to Eq. fl3.1ip , the 



dynamics of the missing bonds is now present for any value of the isospin q. Thus Eq. ( |3.13| ) 
has all the relevant properties of the t- J model we are interested in and we consider it as a 
physically sensible generalization. 
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In the canonical basis the operators ( |3.10|) are given by the following expressions 



d = 2qci4, xu = 4v 2( l- 



Qsi = Cic\(q - a\ai), xl = \2q _ a \ a i c 



Q +i = Cic\\j2q - a\ai a h xu = c\a t , 



Q-i = Cic\a\\j2q- a]a h xl = <4 c i- ( 3 - 14 ) 



Expanding in powers of 1/q the previous realizations, the large-spin limit of Eq. ([3.13 ) turns 



in an interacting spinless fermion hamiltonian imbedded in a classical spin background 
H$ >v] - H* = -2qt^2^(^b^cl Cj + 



'j 



ij 

(3.15) 

where the overlap between even coherent states is 

(n u \n bj ) = (cos | cos I + c -«*-*> sin | sin |) * , (3.16) 

(so that the factor (f^lf^) 1 ^ 2 ^ in Eq. (|3.15|) is g-independent), and we have chosen the 
relative normalization of the multiplets (2.5) such that u>i — (pi/2 = 0. 

A relevant feature of the effective hamiltonian Eq. ( 3.15|) is that it gives a variational esti- 



mate of the ground-state energy. Any further inclusion of 1/ q- fluctuations can only improve 
the estimate, and this is a remarkable property that in our approach is naturally preserved. 
To show this property, we consider the model Q3.13 ) in the original frame, namely with the 



operators (P^Aq) directly replaced by (Q^,Xa), and denote it with H^ Q y According to 
Eqs. 02.151 , 12.161 ) and to the generalization ( |3.3| ) of the charge and spin operators ( |2.2b| ) the 
action of the hamiltonian H^q, over the state 

N h / N \ 

I* (9) > = I1 E^kl |n>, (3.i7) 
i=i \i=i ) 

is exactly given by the action of the effective spinless fermion hamiltonian ( |3.15| ) over the free 
particle state obtained by replacing in Eq. (|3.171 ) the operators xl with spinless fermion 
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operators c] and the boson state with the vacuum state |0) of the spinless fermion 
representation. Here |f2) = ® \^u) is the tensor product of even coherent states and ipi{i) 
are Nh complex and orthonormal orbitals (I = 1,2, ...,Nh in a d- dimensional hypercubic 
lattice L d = N). 

As a consequence, any variational estimate £ = ($\H£j\ty) of the ground-state energy 
of the spinless fermion hamiltonian ( |3.15| ) is also a variational estimate of the ground-state 
energy of the generalized t-J model ( |3.13j ). Moreover simple algebra shows the notable 
identity 

£ = (¥^\H^ Q] \¥^) = 2q (¥^\H tJ \¥^), (3.18) 

where H t j is the exact q = | model ( |2.1|) . The factor 2q simply rescales the unit of energy and 
this ensures that no spurious phase transitions as a function of q are present in connecting 
the large-spin "classical" solution £ to the variational estimate of the exact model. For the 
time being we thus set 2q — 1. 



IV. LARGE-SPIN MEAN-FIELD 

Despite the notable simplification given by the disappearance of the constraint of double 
occupancy, Eq. (|3.15|) is still a difficult problem to solve because the spinless fermions interact 
via a four-body term. For the infinite-f/ Hubbard model, i.e., when J = 0, and for the case 
of a single hole the interaction term does not play any role. Thus in the large-g limit we end 
up with the study of a free hamiltonian with a complicated site dependent hopping, which 
can be easily solved numerically on a finite - but large - lattice, as we shall discuss. For 
low-doping the Hartree-Fock factorization ( |3.17| ) is expected to be a good approximation 



since few holes can rarely interact. We extended the Hartree-Fock solution to the full phase 
diagram of the hamiltonian ( |3.15| ). Although this is not completely justified the interaction 



between spinless fermions should not drastically affect the mean-field phase diagram, since 
for example it is known that the the Hartree-Fock solution is exact in d — > oo for a gas of 
interacting spinless fermions 0. 
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Thanks to Eq. (|3.18|) it is possible to apply Wick's theorem and the classical energy on 
a given arbitrary state reads 

£ = Qij[—t(Qbj\^bi) 9i 



>',j 



- + 9jj - 9i,i9jj + 9i,j9j,i ~ 1)]; 



(4.1) 

where = E£i ^f(«M0')- 

The classical energy is therefore a function of the 2 A" spin angles 9i, (pi and of the 
Nh x A" complex variables defining the wavefunction Q3-17 ). The lowest possible energy 



So as a function of all the 2A" x (1 + Nh) real variables represents a variational classical 
estimate of the ground-state energy. We solved this optimization problem on a square lattice 
numerically, following a scheme which is similar to the one introduced by Car and Parrinello 
[T9| for the electronic structure problem, i.e., for the simulation of the slow dynamics of 



the ions interacting via a self-consistent potential generated by the electronic fast degrees 
of freedom. This approach can be extended to our case because the isospin q behaves as an 
adiabatic parameter: in the q — > oo limit the spl(2,l) even generators ( |2.7| ) ("slow variables") 
become classical objects, so that Eq. ( |3.13p describes the dynamics of the odd and projection 



operators ( 2/7) , (|3.9| ) ("fast variables") in the background of the former. Thus the electronic 



degrees of freedom have a much faster dynamic of the spin angles, which in this case play 
indeed the role of the ionic coordinates. 

In order to minimize the numerical effort we first move the electronic degrees of freedom 
at a fixed spin configuration, and then the spin angles without changing the electronic 
degrees of freedom. At each step we do not require a fully converged Hartree-Fock solution 
of the electronic part, but we make a fixed number of steepest descent steps, followed by 
a Graham-Schmidt orthogonalization of the orbitals to achieve numerical stability. After 
many spin-moves (~ 10, 000) we get a fully self consistent solution of the problem up to 
computer machine accuracy. The solution found with the previously described iterative 
scheme may not coincide with the absolute minimum - the classical ground-state. The true 
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minimum can be identified on a reasonable ground by resorting to symmetry considerations 
or by performing several simulations with different random initializations. 



For J = the mean-field hamiltonian ( |3.15| ) becomes the free spinless fermion hopping 
hamiltonian 

H f = -t^iA^bj)^ (4-2) 

ij 

which has been phenomenologically introduced |lTJ to study the instability of the Nagaoka 



state - i.e., the free particle state with maximum allowed spin and lowest kinetic energy 
- in the infinite-f/ Hubbard model. The Nagaoka theorem |2(J states that it is the exact 



ground-state of the model in the one hole case. Writing the coherent state overlap in the 



form [17 1 



1 + nfij aa 



(fl K |n 6i ) = Y y J e iA " (4.3) 

where Ajj is the solid angle subtended by the vectors Hi, Hj, and a reference fixed one Ho, 
it has been conjectured that for more than one hole a non-fully polarized spin background 
may have better energy than the Nagaoka state, because the gain in magnetic energy due 
to the nonzero flux of Aij could overwhelm the narrowing of the effective bandwidth tQij — > 



tQijyJ (1 + Hifij)/2. In this case we have performed over 10, 000 fully converged minimizations 
on a 10 x 10 and a 16 x 16 lattice and for arbitrary density. We initialize randomly the spins 
angles and set the orbitals at the corresponding Hartree-Fock solution of the spinless- 
fermion hamiltonian with the chosen random spin configuration. In all the minimizations we 
have never found a solution with corresponding energy lower than the "Nagaoka energy". 
Most of the runs converge to the fully polarized solution. Only a few remain trapped 
into a local minimum of the classical energy. A particularly interesting case is at doping 
5 = Nh/N = 1/2. In this case we have indeed found a stable planar solution with fa = 
and Qi = ^r% x , where i x = 0, 1, ...,L — 1 are the lattice coordinates along the x-axis. This 
minimum configuration is similar to the one proposed by Doucgt and Wen for few holes [|TT] . 



However we have found that this kind of state is always unstable except for this particular 



doping 5 = 1/2, and its energy is only degenerate with the Nagaoka energy. For large L 
this state tends (locally) to the Nagaoka state and leads in fact to the same correlation 
functions. Contrary to the Nagaoka state, this planar solution is not an exact eigenstate of 
the infinite-?/ Hubbard model. This fact represents a very simple proof that the Nagaoka 
state is not the true ground-state at 5 = 1/2. Although our result about the stability of 
the Nagaoka state is based on a numerical optimization problem which is never completely 
reliable, we at least may argue that there is no evidence of any exotic spiral, flux or chiral 
phase in the t-J model for J = at the mean- field level. 

For J ^ we consider first the interesting case when only one hole is present. A recent 
reason of debate is whether the spiral mean-field solution of Shraiman and Siggia has lower 
or higher energy compared to the polaronic solutions (see Fig. p]). As we have anticipated, 
also in this case we are left with the study of a free hamiltonian with a complicated site- 
dependent hopping, to be adjusted as to minimize Eq. ( |4.1| ). The results of our numerical 
investigation confirm that the polarons are the lowest-energy configurations of the large-spin 
limit of the t-J model. The small-polaron solution, also called five-sites polaron (the hole is 
trapped in a given site and its nearest neighbors: see Fig. [TJa), is stable for J/t > 0.243. For 
smaller J/t the size £ of the polaron is gradually increased (the 8-sites polaron in Fig. |T]b 
is stable for 0.148 < J/t < 0.243) until £ ~ J~K consistent with theoretical arguments ||, 
so that in the J/t — > the polaron solution eventually turns in the Nagaoka state. With 
the polaron spin background the mean-field energy and the orbital can be easily evaluated 
analytically. We report the case of the five-sites polaron localized at site i . Writing the 
orbital in Eq. ( |3.17| ) in the form = f5i i + giQ ioi , the diagonalization of Eq. Q3.15| ) gives 

£ = -^(2iV-4)-2iA , 



Ao = -¥ + sign(t) V 1+ V8ty 

/ = ~A=s 9, = (4-4) 



The energy of the antiferromagnet with four missing bonds is recovered in the J/t — > oo 
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limit, because we have Ao — > so that / — > 1, gi — > 0, which means that the hole is 
statically placed on site %q. The one-hole energy (referred to the half-filled ground-state 
energy £af = —JN) for the Shraiman and Siggia state at momentum k = (f , f ) and the 
polaron estimate is plotted in Fig. ^| as a function of J/t in a 10 x 10 lattice. The comparison 
of the two energies is meaningful because both represent a variational estimate of the ground- 
state energy of the t- J model. However the two variational wavefunctions refer to different 
order of approximation, because the hole in the Shraiman and Siggia estimate has a definite 
momentum. Propagating the polaron through the whole lattice via processes as depicted in 



Fig. y (highly non-perturbative in 1/q, as discussed in Ref. 0]), eventually decreases the 
energy, but the present data shown in the Fig. clearly indicate that the spiral has lower 
energy than the static polaron, unless for small J/t. 

We consider next the case of four holes. For large J/t the stable state is clearly localized 
in a 2 x 2 plaquette (see Fig. |I]c). As we decrease J, for J/t < 1.09 the four holes prefer to 
remain localized in a polaron of 13 sites, tilted by 45 degrees respect to the x, y axes (see 
Fig. |I|d), rather than to split into a couple of bound pairs. For J/t smaller and smaller the 
size of the polaron smoothly increases until the spins are fully polarized in the given finite 
10 x 10 square lattice. These results support the Emery et al. scenario. At least at the mean 
field level the classical solution consists of an hole-rich phase fully polarized and a classical 
antiferromagnetic region which are completely phase separated for arbitrary J. The exact 
phase boundary in the thermodynamic limit may be obtained by minimizing the energy of 
the phase separated phase. Following Emery et al. we get a critical J c above which phase 
separation occurs, given by 

Ep 

I (E F ~ E)N(E)dE 
T = ~ 2B ' <«> 

E F 



5 = j N(E)dE, N(E) = ^r^l - (E/Atf ); (4.6) 

-A\t\ 

where the spinless fermion density of states N(E) is expressed in terms of the complete ellip- 
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tic integral of the first kind K, Ep is the spinless fermion Fermi energy at the corresponding 
doping 5, and B is the classical energy per bond (B = 1) of the Heisenberg antiferromagnet. 
The resulting phase diagram in the J — p plane, where p = 1 — 5, is shown in Fig. [|. Here 
we note the characteristic 5 ~ yfj singularity occurring at low doping as the size £ ~ of 
the polaron increases in a singular way. At low density Eq. ( |4.5| ) is valid only if electrons do 
not form bound states. This is actually the case for the exact t- Jhamiltonian ||, where two 
electrons form a bound singlet pair for J/t > 2. In this approach however the paramagnetic 
phase at low density is clearly not well characterized since the spins are frozen in some fixed 
direction. As a result two electrons never bind - as we have tested numerically - and the 
Emery-like phase boundary represents the exact mean-field phase diagram. Let us discuss 
now what happen for J < J c at fixed density. At J = J c the Nagaoka state is stable by 
construction since the hole-rich phase exhausts all the allowed space at the given density. 
It is then easy to convince ourselves that if we decrease J/t the Nagaoka state is even more 
stable because lowering J/t favors the polarized solution. 

At the large-spin mean-field level we therefore end up with a very simple phase diagram 
consisting of a ferromagnetic phase for J < J c and a separated phase for J > J c . 

V. DISCUSSION AND CONCLUSIONS 

In this paper we have presented the coherent states that allow to give a variational 
foundation to the large-spin limit of the t- J model. These states are obtained by grading ac- 
cording to the spl(2,l) algebra the standard su(2) spin coherent states. We have emphasized 
how the large-spin limit can be performed in distinct schemes, and applying one of them 
we have rigorously derived the Kane, Lee, and Read hamiltonian. We have discussed the 
shortcoming of this approach, related to the fact that in this scheme the mean-field solution 
turns out to be the classical Neel state, independently from the density. As a consequence 
phase separation cannot be recovered at the mean-field level. 

We then have presented and numerically investigated the effective spinless fermion hamil- 
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tonian obtained by means of a different and more satisfactory definition of the large-spin 
generalized t- J model. This mean-field model hamiltonian has the essential properties and 
symmetries of the original hamiltonian, in particular it exactly takes into account the con- 
straint of no doubly occupancy. The mathematical tools we have employed are similar to 
other proposed techniques, notably the Swinger bosons/slave-fermion representation. How- 
ever we hope to have better clarified the mathematical and physical assumptions underlying 
the large-spin limit of the t-J model. More important, by means of the spl(2,l) graded 
coherent states we have shown that the mean-field solution is independent of the magnitude 
of the expansion parameter and therefore the corresponding energy is variational. In this 
way we have generalized a very useful property of the spin systems. 

Our numerical investigation confirms that the polaron solutions are the lowest-energy 
configurations of the t-J model in the large-spin limit. However we do not find that the 
static small-polaron has lower energy compared to the Shraiman and Siggia spiral state 
as reported by Auerbach and Larson ||21|| , at least when in the definition of the t-J model 
three-site contributions are neglected. The investigation of the full phase diagram at the 
mean-field level support the Emery et al. picture: in the J — 5 plane a critical J c = J c (5) 
with J c — > for 5 — > separates the ferromagnetic and the phase-separated phases. 

For J = the large-spin solution is always fully polarized and our numerical analysis 
suggests that the instability of the Nagaoka state cannot be found at the mean-field level. 
The mean-field picture is exact for a single hole, likely for small doping, but is incorrect 
for large doping, where existing numerical works and the variational singlet Gutzwiller 



projected wavefunction [^j result seems to suggest that the Nagaoka state never survives 
at finite concentration of holes. On the other hand for J = the energy So obtained by 
minimizing the functional ( |4.1| ) has the particle-hole symmetry, i.e., Sq(5) = £o(l — This 
is not a true symmetry of the infinite-t/ Hubbard model and thus we expect that 1/q- 
fluctuations play a relevant role, especially for small density. Non-perturbative processes 
may be of some relevance, however we expect the mean-field picture to be reliable in the 
large J/t region. 
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FIGURES 

FIG. 1. Polaronic solutions. The symbols o and • indicate spin up and spin down, respectively. 
The symbols <S> and indicate spin up and spin down, and in addition a nonzero hole density on 
the site. For the one hole case we report the small-polaron (a), the 8-sites polaron (b), and the 
13-sites polaron (d). For the four holes case it is reported the 2x2 plaquette (c) and the 13-sites 
polaron (d). 

FIG. 2. Variational hole energy in unit of \t\ as a function of J on a 10x10 lattice. The 
dashed line is the semiclassical spiral solution of Shraiman and Siggia with momentum k = (§,§)• 
The continuous line represents our large-g mean-field solution: a static polaron of increasing size 
as J is lowered. The abrupt changes of the energy at small J are due to corresponding first order 
transitions to larger and larger polaron size. 

FIG. 3. Second order process for the q = \ theory that allows the small-polaron propagation. 
The wavefunction of the polaron at site % is written as = \ fi) + J2k ^ik\9ik)- I n ( a ) we nave 
the state |<7j&) entering the polaron at site i. In (b) we have the intermediate state \5j) ~ H t j\gik) 
which is a one-particle excited state and corresponds to a local minimum of the energy functional 
(4.1). In (c) we have the state \fj) entering the polaron at site j, which is obtained by a further 
application of the hamiltonian: \fj) ~ H t j\Sj). Due to (^i\H^j\^j) ~ giftJ, for % / j, an effective 
intrasublattice hopping is generated. 

FIG. 4. Phase diagram of the t- J model obtained with the Hartree-Fock large-spin approach. 



26 



